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ДИНАМКА КОРЕНЕВОГО ПОРТРЕТУ 1НТЕРВАЛЬНО! 
СИСТЕМИ НА МЕЖТ СТИКОСТИ 


В статье рассматривается задача определения характера распределения значений параметра 
корневого годографа интервальной динамической системы вдоль границы устойчивости. 
Используются поля корневых траекторий. В результате обеспечивается возможность установления 
значений параметров, при которых система сохраняет устойчивость. 
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У стати розглядаеться задача визначення характеру розподмлу значень параметру кореневого 
годографа 1нтервально! динам!чно! системи вздовж кордону стЙкост!. Використовуються поля 
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при яких система збертгае стйюсть. 
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Существенное значение при синтезе качественных,  робастных 
динамических систем имеет анализ действующих на систему возмущений и 
вызываемых ими неопределенностей [1—6], установление характера влияния 
неопределенностей на динамику системы. 

Рассмотрим динамическую систему с неопределенностью, описываемую 
семейством интервальных [3—6] характеристических уравнений (полиномов) вида 


п . 
п-) _ 
Ха, °=0, (1) 
1=0 
где вектор коэффициентов а = (ао, ар, ..., ав, а.) принадлежит некоторому 


связному множеству А < В", а е А, п — степень полинома (целое число); $— 
комплексное переменное, $ = с + 1. 
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Интервальное изменение коэффициентов — характеристического 
полинома системы отражает интервальное изменение ее параметров. 
Предположим, что коэффициенты (Т) действительны и изменяются в 
следующих пределах: 


а; <а, <ау, ]=0,п, ау = 0, (2) 
где а; и а; — соответственно минимальное и максимальное значения замкнутого 


интервала (2) изменения коэффициента а. 

Одним из эффективных подходов к исследованию данного вопроса является 
корневой подход [5-9], который по своей природе предполагает вариации 
параметров и является мощным средством синтеза и анализа, позволяющим 
представить полную картину изменений, которые происходят в системе в ответ на 
параметрические вариации. 

Сделав в выражении (1) замену переменного, 5=& +1, запишем его 
относительно некоторого коэффициента а», который назовем варьируемым 
коэффициентом: 


$($) 
а; = 7(5) Ро Е и(д, @)+1’(о, @)), (3) 
(5) 
где и и у - гармонические функции двух независимых действительных 
переменных оси @; $5) и $) - некоторые полиномы от комплексного 


переменного 5. 

Выражение (3) представляет собой рациональную функцию отображения 
интервального семейства подобных функций, которая позволяет отображать 
некоторые образы, заданные в плоскости варьируемого коэффициента а», на 
плоскость комплексного переменного 5, т.е. фактически является уравнением 
корневого годографа системы в общем виде. 


Запишем функцию параметра годографа [5, 7], 
а; =и(0, 0), (4) 
и уравнение корневого годографа Теодорчика-—Эванса [5, 7], 
у(о, 0) =0. (5) 


Исследуем поведение функции параметра корневого годографа (4) на 
границе устойчивости — оси 1% плоскости корней 5. С этой целью рассмотрим 
семейство Р двухпараметрических полей, К; е РЁ, свободных корневых траекторий 
Теодорчика — Эванса [5] системы с параметром поля а„ 1 и исследование будем 
выполнять для отдельного поля А; (1 = 1, ..., 0) как характерного представителя 
интервального семейства полей, при условии —0< а„ < +50 ‚ —0< а, 1 < +00. 


Вначале, подставив значение о’= 0 в уравнение (4), определим для поля А; 
распределение значений функции параметра вдоль оси 14 плоскости корней: 


— при и четном, 7/2 четном 
и( ©) =— а@" + а@"? — аао"\ +... але‘ На @?=а,; (6) 
— при и четном, и/2 нечетном 


и( ©) = ао@ — а2@"? + ад“ -... - ар дай 4 а,2@? = ан; (7) 
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— при и нечетном, (и — 1)/2 четном 
и( 6) =— а! + аз" 3 — а@-... - ав 4 + ара? = ая; (8) 
— при и нечетном, (и — 1)/2 нечетном 
и( ©) = а"! — аз" 3 + а; 5-... а 4 + а.о@? =а,. (9) 


После подстановки значения о’ = Ов (5) определим уравнение корневого 
годографа для оси №: 


— при п нечетном, (и — 1)/2 четном 
Хо) = а — а@? + а." 4... - аз +а,л@= 0; (10) 
— при и нечетном, (и — 1)/2 нечетном 
4) =-— а0@ + а2@" ? — ал 4+... - аз +а@= 0; (11) 
— при п четном, 7/2 четном 
хо) =-а@ 1 + аз 3 — +... - аз +ае= 0; (12) 
— при п четном, 7/2 нечетном 
Хо) =а1@' ! — аз" 3 + аз °-... аз +а,ло= од. (13) 


Для определения значений координат @ в точках пересечения границы 
устойчивости ветвями корневых годографов поля используются уравнения 
корневого годографа (10) — (13). Для определения значений свободного параметра 
(а) в найденных точках применяются выражения (6) - (9). 

По причине свойств аналитичности и непрерывности функций (6)- (13) 
очевидно, что при пересечении годографами поля корневого портрета границы 


устойчивости (оси 16) на этой оси формируется некоторая многолистная область 
БУ, и каждая ветвь каждого годографа образует непрерывную подобласть (лист) 


этой области. 

Определение 1. Областью пересечений поля корневых траекторий назовем 
область пересечений границы асимптотической устойчивости динамической 
системы ветвями корневых годографов поля корневых траекторий этой системы. 

Заданный интервал изменения параметра поля а„,! определяет значения 
границ реальной области пересечений поля. Для полного поля (- © < а, |< +00) 
область пересечений охватывает всю границу устойчивости. Множество 


координат 7, области пересечений определяется как множество И. 
Ы. Л 
@р, Е ЙЪ, 


где ©), — значения координат точек области р/ пересечений границы 


устойчивости ветвями корневых годографов поля; Й/ — множество координат 
Л Л Я тд Л 
@р, точек области О.,. Множество значений параметра ар, в точках @р, 


обозначается как А? са 


Я Е рр 
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С целью определения максимального и минимального значений функции 
параметра в области БУ выполняется ее исследование на экстремум при помощи 


первой производной. Производная 1’(®) от функции параметра (6)- (9) 
приравнивается нулю: 


п четное, и/2 четное: 
и’(@) = паи" 1 (и-2)а2 "3+ (п 4)ал а 5-...+ 4а, аа -— 2ао@=0. (14) 
п четное, и/2 нечетное: 
и’(6) = па" (п-2)а› "3 + (п-4) ал" 5-...— Ча, аа + 2а,о®=0. — (15) 
п нечетное, (п — 1)/2 четное: 
и’(®) =— (и-Па "2+ (п-З)аз@" — (п-5) аз" 6+.. — Ча, ла + 2а,о@=0. (16) 
п нечетное, (п — 1)/2 нечетное: 
и’(6) = (и-Гаа?- (п-З)аза" “+ (и-5) аз" 6-...— 4а, а 2а,о®=0. (17) 


Из четырех последних уравнений определяются координаты @/. точек 


экстремума функции параметра на границе устойчивости: 


т /. (18) 


ет е 


где И’/ — множество координат точек экстремума @/. 


Для определения множества А/, экстремальных значений а/, параметра 


в точках (18) используется выражение 


С целью получения экстремальных значений свободного параметра 
Л 
множества и используются выражения (6)- (9), в которые подставляются 


предварительно вычисленные из уравнений (14) — (17) соответствующие значения 


Л 
>; й 
Р. хх 
Максимальное значение параметра а Ри на границе устоичивости определяется 
выражением 
Л = Л 
Я о тах = р 4», (19) 
минимальное значение Параметра вы на границе устойчивости выражением 
Л а Л 
По = ШЕЯ а (20) 
1%. Л Л 
Значение © в точке Ч тах обозначается как > тах › а значение © в точке Чит как 
©’ 


е шш ° 


На рис. 1 представлен пример (эскиз) графика распределения функции 
параметра вдоль границы устойчивости при рассмотрении полного поля [5] 
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корневых траекторий системы. На графике с целью упрощения соединяющие 
вершины кривые изображены прямыми линиями. Знаком «+» на рис. 1 
обозначаются участки возрастания функции, знаком «—- убывания. В 
соответствии со свойствами корневых портретов систем, приведенными в [5], при 
положительных значениях коэффициентов характеристического уравнения 
системы на графике, начиная от нулевой точки по направлениям вверх и вниз, 
вначале располагаются участки возрастания (при значении параметра поля а, 1 =0 
в начале координат располагаются две начальные точки, на правой 
действительной полуоси начальные точки отсутствуют). Далее имеет место 
чередование знаков соответствующих участков. Чередование участков имеет 


место в пределах определенной области Ё у на границе устойчивости, называемой 


областью экстремумов. За пределами области экстремумов функция имеет 
монотонно возрастающий/убывающий характер. 

Определение 2. Областью экстремумов функции параметра корневого 
годографа на границе устойчивости системы назовем область, в пределах которой 
располагаются точки экстремума (19) - (20) этой функции. 

В пределах поля корневых траекторий область экстремумов определяется 
максимальной 6 . 


и минимальной 6, границами (рис. 1), которые вычисляются 


из следующих выражений: 


©/ = зар И/ 


@/ = ШЕИ. 


Таким образом, область пересечений р; при рассмотрении полного поля 
разделяется на три следующие области (рис. 1): 


— область экстремумов Ё 7 ‚ определяемую в пределах 


— область возрастания/убывания функции параметра верхнюю 7 ыы 


определяемую в пределах 


@) <) <+%; (21) 


е 


— область  возрастания/убывания функции параметра нижнюю р: Сы 


определяемую в пределах 
а: (22) 


Следует отметить, что реальные границы области пересечения границы 
устойчивости ветвями годографов заданного поля, определяемые интервалом 
вариации параметра поля, и границы этой области, определяемые выражениями 
(21) - (22) для полного поля, могут не совпадать. 


© А.А. Несенчук 55 


155М 1561 - 5359. Штучний 1нтелект, 2015, № 1-2 


10 


Ц > 
Че алии я 
5 
7 
Че; тах 
+ ` 
у 3 
Я =@ 2 [а] 
== + 2, => 
-4 с 
# 
Яо. 
Е: 
Че: 
> © 
я 


Рис. 1. Распределение значений функции параметра вдоль границы устойчивости. 
Область экстремумов 


При задании определенных границ интервалов вариации параметра поля 
(когда поле не является полным) область пересечений р/ не бесконечна, т.е. 


—. . 


является ограниченной. Максимальная › И минимальная ©@ъ координаты, 


ограничивающие область пересечений, определяются из уравнения корневого 
годографа (10) - (13), которое описывает динамику координат точек пересечения 
на границе устойчивости. 


Перепишем уравнение корневого годографа, например (13), относительно 
параметра поля аи: 


— во"! + ад" 3 — 4+... +а,за? = а. (23) 


Функция (23) является аналитической непрерывной аналогично функциям (6) — 
(9), и характер распределения значений параметра поля а„1 вдоль границы 
асимптотической устойчивости аналогичен характеру распределения значений 
свободного параметра а„, которое было рассмотрено выше (рис. 1), т.е. представляет 
собой чередование возрастающих и убывающих участков на соответствующем 


графике. Следовательно, граничные значения 4), © и © ‚ соответствуют 


предельным значениям интервала вариации параметра поля а„.1, и для определения 
реальных границ области пересечений ветвями годографов поля границы 


устойчивости необходимо найти конечное множество И! ‚ координат этих 


граничных значений (для каждого участка области пересечений в общем случае их 
будет два). Для этого необходимо решить уравнение (23), подставив в него 


предварительно вначале максимальное @а,„_|: 
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— ао" |+ да 3 — 45+... на: =а,, 
а затем минимальное аи-—1: 
— ао"! + а" 3 — 5+... а за= аи 


значения заданных интервалов изменения параметра Поля и выбрав из числа 


полученных корней ©, максимальное @/ и минимальное © значения: 
= Г 
р == 5ир . 2’ 


@1 = ШЁ И. 


Аналогично следует поступать и в случаях использования уравнений (10) — (12). 


Рис. 2. Область пересечений границы устойчивости 1@ ветвями годографов поля 
корневых траекторий системы шестого порядка. Реальная область пересечений 


Пример реальной области пересечений, лежащей в пределах области 
экстремумов, приведен на рис.2, где представлено поле корневых траекторий 
динамической системы шестого порядка. Полная область пересечений при 
изменении параметра поля а„ 1=а5 в пределах 0 < а5 < +0 на этом рисунке обозначена 


как 7 . Реальная область пересечений, формируемая при 0<а5<а5 (где а5-— 
некоторое заданное значение параметра поля), которую обозначим как ру Е 
(имеется в виду ее часть, расположенная в верхней полуплоскости) состоит из 
подобластей Б,. . РУ и В формируемых при пересечении границы 
устойчивости одной из двух пересекающих эту границу ветвей каждого годографа 


поля, т.е. ру = И , +07 5+... Годограф, показанный штриховой линией, не 


рассматривается, поскольку он полностью не устойчив. 
Следует отметить, что границы © и © реальной области пересечений 
заданного поля могут различным образом располагаться по отношению к 
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границам области экстремумов. Например, реальная область пересечений может 
включать в себя всю область экстремумов или только ее некоторую часть (рис. 2), 
а также располагаться за пределами области экстремумов, т.е. быть частью 


рассмотренных выше областей р/ гой РУ,, лежащих за пределами области 


@Ь› 
пересечений. 

Буквами 41—44 отмечены граничные точки диаграммы (заштрихованные 
участки) распределения значений свободного параметра а, годографа вдоль 
границы устойчивости в пределах реальной области пересечений для заданных 
границ изменения параметра поля, а, 1 = 45. 

Для определения экстремальных значений функции параметра корневого 
годографа на границе устойчивости для всего интервального семейства 
годографов корневого портрета системы необходимо в выражения (6)- (9) 
подставить соответствующие значения границ интервалов изменения 
коэффициентов: от а1 до а» 2. 
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ВЕ5ОМЕ 
А.А. МезепсвиК 
Оупапис$ оГ Ме Пиегуа! Сопиго! Зузет Вооё Госи$ Роггай аё Фе За Иу 
Воипдагу 

Тре рарег 4еа]5 ул шуезИеайоп оЁ фе ргоет оЁ ащютайс сопёо| зузеп$ 
тоби$ез$, 1.е. йе еНесЕ оЁ Фе р|апЕ рагатеег$ забзапна| аеуланоп$ (уапаНопз) ут ве 
сопёго! зузет згасвте, саизше зузет рагатейлс ипсенайиу, ироп фе зузет 4упапис$, 
15 а БИ ап ореганоп аиаШу. Рупапис зузетл$ у пщегуа| рагатебле ипсенаштиез 
Вауе Бееп сопз14егед. шуезисаноп оЁ фе зузетаз Бебауюг ш соп 01$ оЁ рагатейлс 
уапанопз Ваз Бееп сагле4 оп. Рог 15 ригрозе Ше ТеодогсК — Е\’апз Нее гооф 10си$ 
ройгайз апа Не!45 Вауе Бееп аррПед. ТВе зуует гос 1осиз ропгай 15 гергезещеа Бу Фе 
ЁатйПу (Фе $60) оРбмо-рагатеЕг гоо! 1оси$ Не!4$ шБепипе зпип|аг ргорегиез \асВ дебпе Фе 
ргорегйез оЁ Ве зу$ет аз а \пое. Тре рипстр[ез ог сепегайпе йе Фаотат оЁ Фе зубет 
гоо{ 1осиз рагате{ег ФЧзитиНоп аоп» Фе за БИ Боппдагу Юг Фе гос! 1осиз Йе]4 вВауе 
Бееп оМашед, усе аге аррИе4 юг Нпаше пиегуа[$ оЁ рагатеег$ епзипие фе зузет 
гоби$Е заб у. 

Рагфег суа[ланоп оЁР фе уогк со Ъе несе ю \огкте оп Ше 1есмаие Юг 
са1сШайпте 4упапис зузет сБагасепзис едианоп соеЁйсет® ап4 рагатеег$ епзайие 1$ 
тоби$Е заб ИШу ап гоби$Е диаШу, Базед оп фе 41зсоуеге4 теошапиез оЁ Ше уапаЫе зузет 
рагатеетз уааез Язитбаноп а1опо йе за у Боипдагу. 


А.А. Несенчук 
Динамика корневого портрета интервальной системы на границе 
устойчивости 

Статья посвящена изучению вопроса робастности автоматических систем, а 
именно, влияния существенных отклонений (вариаций) параметров объектов 
управления в структуре систем автоматического управления, вызывающих парамет- 
рическую неопределенность, на динамику системы, ее устойчивость и качество функ- 
ционирования. Рассмотрены динамические системы с параметрической неопре- 
деленностью интервального характера. Выполнено исследование поведения систем в 
условиях параметрических вариаций. С этой целью использованы свободные корне- 
вые портреты и поля корневых траекторий Теодорчика — Эванса. Корневой портрет 
системы представлен семейством (множеством) двухпараметрических полей корневых 
траекторий, обладающих одинаковыми свойствами, определяющими свойства всего 
портрета в целом. Определены принципы построения диаграммы распределения 
значений параметра корневого годографа для поля корневых траекторий динами- 
ческой системы вдоль границы устойчивости. На основании подобной диаграммы 
могут быть установлены интервалы изменения значений параметров системы (интер- 
валы неопределенности), при которых данная система является робастно устойчивой. 

Дальнейшее развитие работы может проходить в направлении разработки 
методики расчета коэффициентов характеристического уравнения и параметров 
динамической системы, обеспечивающих ее робастную устойчивость и робастное 
качество, основанной на установленных закономерностях распределения значений 
варьируемых параметров системы вдоль границы устойчивости. 
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